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ABSTRACT 

Using new and published photometric observations of /j, 1 Sco (HR 6247), spanning 
70 years, a period of 1.4462700(5) days was determined. It was found that the epoch 
of primary minimum suggested by Shobbrook at HJD 2449534.178 requires an adjust- 
ment to HJD 2449534.17700(9) to align all the available photometric datasets. Using 
the resulting combined-data light-curve, radial velocities derived from IUE data and 
the modelling software Phoebe, a new system solution for this binary was obtained. 
It appears that the secondary is close to, or just filling, its Roche-lobe. 
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1 INTRODUCTION 

fj, 1 Sco (HR 6247; HD 151890; HIP 82514) was only the 
third spectr oscopic binary to b e discovered and is listed in 
the GCVS (|Samus et alj|2003 ) as an eclipsing binary vari- 
able. Over the intervening years there have been a number 
of detailed measurements and several significant studies. It 
is now classified as a semi-detached (sd) binary as one com - 
ponent is believed to fill its Roche lobe. ICester et al.l (|l977l ) 
considered fi 1 Sco, to be unusual owing to: 

(i) their determination of an apparently high mass-ratio; 

(ii) indications that the secondary component has a larger 
radius than the primary and overflows its Roche lobe; 

(iii) both components appear to lie on the main sequence; 

(iv) it being a member of the Scorpius-Centaurus cluster 
thus indicating an age of no more than 10 7 years hence the 
possibility this system has just arrived on the main sequence. 

A comprehensive analysis of fj} Sco was undertaken to 
determine the best estimate of period from all available pho- 
tometric data, establish a suitable epoch and calculate the 
key parameters defining the system. 

The photometric data used included published measure- 
ments from: 



(i) IRudnick and Elyevl i|l938h . 
(if) Ivan Gentl dl939l ), 

(iii) IStibbsl l| 19481 ). 

(iv) Tycho collected be tween 1990 and 1993 

^European Space Agencvlll997l ; lOchsenbein et alJlioobl ) , 

(v) Hipparcos collected between 1990 and 1993 
^European Space Agencvlll997l ; lOchsenbein et alJlioobT ), 



(vi) IShobbrookl (|2004l ). 

and new measurements taken in 2006, 2007 and 2008 by one 
of the authors (Moon). 

The photometric data was then combi n ed w ith ra- 
dial velocities obt ained by Stickland et al. (Il996f) from 
IUE spectral data ((international Ultraviolet Explorer Datal 
2009) and the resulting d ataset analy se d using the soft- 
ware program PHOEB E 1 Prsal 20031 : IPrsa and Zwitterl 

120051 : IPrsa et all 120081; iPrsal 120091) based on the estab- 
lished IWilson and Devinnevl i|l97ll ) theoretical construct 
jKallrath and Milonelll999l ). 



2 PUBLISHED DATA 

The sta r p, Sco was first reported t o be a spectroscopic bi- 
nary bv lBailevI (|l896l ). lMaur^l (|l920h undertook an extensive 
study of this star using 184 spectra collected over the pe- 
riod 1892 to 1918 and obtained a radial velocity curve of one 
component with respect to the other. From these data she 
determined a period of 1.44627 days and selected an epoch 
for the zero-point of the recession of the brighter component 
of the binary at HJD 2412374.434. The total light variation 
was between 0.3 and 0.4 of a magnitude and the combined 
mass of the system was estimated by her to have a minimum 

of 1 6.5 MfT> 

IRudnick and Elvevl (|l938l ) made 128 photoelectric mea- 
surements of jj 1 Sco using a photometer bas ed on a Kunz 
KH photocell as described bv lStebbinsl l|l93ll) . This cell had 
peak sensitivity near 460 nm and was typically used in con- 
junction with blue and yellow filters to define band-passes 



2 C. van Antwerpen and T. Moon 



with effective wavelengths around 400 nm and 500 nm. Rud- 
nick and Elvey do not, however, indicate the band-pass for 
their measurements or whether filters were used. Using their 
photomet ric mea s ureme nts with the spectroscopic elements 
given by iMaurvl l|l920f) . they determined the mass of the 
system to be 24.0 Mq, consistent with the minimum value 
of 16.5 Mq determined earlier by Maury. Additionally, they 
estimated the effective temperatures of the components to 
be 17,000 K and 13,000 K. 

Around the same time the l ight variations of fi 1 Sco 
were also measured by Ivan Gent) (|l939h using photographic 
techniques. His 86 m pg measurements showed variations for 
the primary and secondary eclipses of 0.30 and 0.19 mag- 
nitudes respectively - similar to that obtained by Rudnick 
and Elvey and consistent with the original estimate made 
by Maury. There was, however, a systematic difference in 
the magnitudes made on the two different systems with van 
Gent's being 0.02 magnitudes brighter than those of Rudnick 
and Elvey (although both sets of measurements used /i 2 Sco 
as the comparison star). This difference was attributed to 
the significantly different effective wavelengths at which the 
two sets of measurements were made as there is only a very 
small difference between the colour indices of the variable 
and comparison star. 

For the magnitude versus phase plot of van Gent there 
also appears to be a small shift in phase when it is compared 
to all the other photometric datasets. This shift corresponds 
to approximately 0.02 of the phase for the secondary min- 
imum which equates to 40 minutes. It is not clear whether 
this is a real feature or a measurement error. When com- 
paring the various photometric datasets (see Figure [TJ his is 
the only dataset displaying this discrepancy. An adjustment 
was thus made to all van Gent measurements along simila r 
lines to that discussed in Mo on and van Antwerpenl (|2009t) . 

In his analysis Ivan Genti |l939j) noted the primary and 
secondary minima appeared to be symmetrically and equally 
spaced. By calculating the difference between the observed 
secondary minimum and the halfway point between consec- 
utive primary minima, he showed the values for the eccen- 
tricity and longitude of periastron given by Maury were in- 
consistent with the observed light curve. Consequently he 
suggested the eccentricity is close to zero and considered a 
circular orbit for his subsequent analysis. This was subse- 
que ntly confirm ed through new spectroscopic observations 
fsee lStibbs 1948). Using the observed light curve, along with 
a radii ratio adopted from Maury's estimate of line intensity 
ratios, van Gent determined the mass of the system to be 
25.5 Mq and the effective temperatures of the components 
to b e 20,000 K an d 12,500 K. 

IStibbsl { 19481 ) results were derived from 358 photoelec- 
tric observations made with a gas-filled sodium cell operat- 
ing in a broad spectral passband from 325 nm to 570 nm with 
an effective wavelength of 397 nm. Using these data, along 
with the ra tio of luminosities of th e two components de- 
termined bv lStruve and Elvevl (|l942l ). he calculated orbital 
elements for the /x 1 Sco binary system and co mpa red his re- 
sults w ith those of I Rudnick and Elvevl (1938) and Ivan Gentl 
(1939) showing that there was general agreement between 
them. Stibbs determined the mass of the system to be 
23.2 Mq; fi 1 Sco was thus considered to comprise a primary 
component around 13 to 14 solar masses with a secondary 
around 9 solar masses in a relative orbit with a radius of a 



little over 10 million kilometres. The secondary was consid- 
ered to be filling its Roche lobe and larger than the primary 
by about 10%. The data tabulated were determined by as- 
suming the measurements made were symmetrical about the 
epoch of the minimum under consideration, and then com- 
bined into normals containing ten observations, the normals 
being formed progressively at intervals of five observations. 
This is in effect making the assumption that the orbit is 
circular. Stibbs then reduced his 358 observations to 72 val- 
ues covering half a period. Only the 72 measurements are 
avai lable from th e litera ture. 

ICester et al.l (|l977l ) analysed data for 12 systems clas- 
sified as being semi-detached binaries that presented what 
was, at that time, regarded as an evolutionary paradox be- 
cause the less massive secondary component appeared to be 
in a more advanced stage of evoluti on than the pr imary. For 
their analysis of fj, 1 Sco they used IStibbsl l|l948h data and 
determined masses and radii consistent wi t h ear lier works. 
In their discussion of this star lCester et al] (|l977l ) remarked 
that the radius of the secondary, which appeared to be over- 
flowing its Roche lobe, coupled with the high mass ratio was 
unusual for sd-binaries. Additionally, they regarded its po- 
sition in HR diagram to be unusual as both components ap- 
peared to be lying on the main sequence. Noting that fj} Sco 
was a member of th e Scorpius-Centauru s cluster whose age 
is at most 10 7 years, ICester et alj l|l977h suggested that this 
star may have just arrived on the main sequence. 

ISchneider et all l| 19791 ) also studied fi 1 Sco using Stibbs 
data and determined masses and radii generally consistent 
with earlier works. They evaluated the effective temper- 
atures of the primary and secondary components to be 
21,500 K and 16,200 K but, in contrast to earlier works, 
found the primary component to be about 5% larger than 
the secondary. 

IShobbrookl (|2004h undertook a program of fresh obser- 
vations of those bright eclipsing binaries which had not been 
observed for many years with a view to establishing more re- 
cent epochs and more accurately determining their periods. 
Observations were made in the Stromgren y band and trans- 
formed to V magnitudes. Shobbrook determined a period of 
1.446270(1) day for fi Sco and established a new epoch for 
its primary minimum of HJD 2449534.178 using his 112 new 
measurements, Hipparcos measurements and those taken by 
Stibbs 17,533 days earlier. 

Hipparcos H p and Tycho Vr magnitudes for fi 1 Sco 
were taken from the Hippar cos and Tycho catalogues 
(|European Space Agencvlll997i) made avai lable through the 
VizieR database (|Ochsenbein et al.ll2000t ). The Hipparcos 
data cover the period from 27 January 1990 to 22 Febru- 
ary 1993, and the Tycho data from 10 February 1990 to 
22 February 1993. In the analysis undertaken, only data 
where the value of the quality flag was less than or equal 
to t wo were used. 

IStickland et al.l (|l996l ) noted the paucity of radial ve- 
locity data available for fi 1 Sco and undertook an anal- 
ysis of spectroscopic data available from the IU E satel- 
lite ([international Ultraviolet Explorer Data! 120091 ) ; in all 
19 measurements for each component. They noted the signal 
of the secondary to be less than one third of that of the pri- 
mary making measurements of radial velocity for both com- 
ponents difficult. Based on their analysis, they concluded 
that the orbit may have a small eccentricity of 0.019 ±0.017. 
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Figure 1. Measured magnitudes of fi 1 Sco as function of phase using an epoch of HJD 2449534.1770 and a period of 1.4462700 day. 
The typical error for a photoelectric measurement is no more than ±0.01 magnitude. For photographic measurements (i.e. those of van 
Gent) the typical error is ±0.05 magnitude. 



Additionally, there could be effects arisin g from tidal in- 
fluenc es, gas streams or gravity darkening. IStickland et alJ 
1 19961 ) determined the system to be at an inclination of ap- 
proximately 62°, with a primary component of 8.6 Mq and 
a se condary of 5 . 6 M^ , just filling its Roche lobe 



lArias et all (120051 ) examined jj, 1 Sco spectra for discrete 
UV features and, in the process, obtained radial velocity es- 
timates from the 18 spectra they collected during July 1990. 
While stating their radial vel ocity estimates to be i n very 
good agreement with those of IStickland et all (|l996i ). they 
only provide a combined spectral plot the resolution of which 
is insufficient for a satisfactory estimate of velocity values 
to be made. Their r esults do, however, confirm those of 
IStickland etafl (,19961 ). 

Figure [5] displays the radial velocity data and resulting 
theoretical fit using PHOEBE. 



3 UNPUBLISHED DATA 

Photoelectric measurements of /i 1 Sco were made by one of 
the authors (Moon) from 3 July 2006 to 8 October 2008 
using an Optec SSP-5A photometer attached to a perma- 
nently mounted 10-cm telescope that is housed in an ob- 
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Figure 2. Measured radial velocities of the components fi 1 Sco 
as function of phase using an epoch of HJD 2449534.1770 and a 
period of 1.446270 day. 



servatory with a roll-off roof. For each measurement, the 
integration time is 10 seconds with each observation being 
the mean of 5 consecutive measurements. As the observa- 
tory is situated in an outer suburb of a major (Australian) 
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Table 1. The 133 new photometric measurements of Moon transformed to Johnson V magnitudes and B — V colour indices. 
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city, the background sky was measured for each observa- 
tion. When measuring through bot h B and V filters the se - 
quencewas Vkar, Bstar, B s ky, Vsky l|Otero and Moonl l2006'). 
The resulting 133 new measurements of V magnitude and 
68 of B — V were determined relative to fj, 2 Sco (HR 6252; 
HD 151985; HIP 82545), the same comparison star used by 
previous observers. The check star, HR 6214 (HD 150742; 
HIP 81972), was the same as that used by Shobbrook. The 
magnitudes and colour indices obtained were transformed 
to the Johnson V band using well-determined transforma- 
tion coefficients. These new measurements are presented in 
Tabled] 

The photometric data on /i 1 Sco used in this analy- 



sis thus represent a comprehensive, and perhaps exhaustive, 
dataset comprising 512 measurements spanning 70 years. 



4 ANALYSIS 

4.1 Correction of photometric measurements to 
Johnson V magnitudes. 

The fj, 1 Sco photometric data assembled by the authors com- 
prises measurements made with different photometric sys- 
tems and is thus a heterogeneous dataset; the particulars 
of the various photometric systems are given in Table [2] 
Corrections were applied to the individual datasets to ad- 
just them to the Johnson V band. Owing to its line-of-sight 
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Table 2. Particulars for the photometric bands of the various datasets used. 



Photometric Band 


Effective Wavelength (nm) 


FWHM (nm) 


References 


Yerkes (KH cell) 


? 460 
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Stebbins (1931) 


Photovisual (m pv ) 
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van Gent C 1939) 
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Allen (1973) 


Johnson B 
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Allen fl973~> 
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397 


120 
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alignment with y} Sco, and similarity in both magnitude 
and spectral type, fj, 2 Sco is an ideal comparison star and 
has been used as such by all observers listed here with most 
fi 1 Sco measurements given simply as a magnitud e difference 
relativ e to fi 2 Sco. In measurements made by IShobbrookl 
and those by Moon, fj, 2 Sco was used as the primary 
comparison star with HR 6214 used as a check star. The 
consistent use of Sco as the comparison star makes the 
adjustment of the various datasets to the Johnson V band 

easier. 

The GCPD (|Mermilliod et all 1 19971 ) value of V=3.565 
was chosen as the magnitude for fj, 2 Sco to which the dif- 
ferences measured for (j, 1 Sco were then applied. Colour cor- 
rections were not applied as there is only a small difference 
0.002 magnitude in the B — V colour indices of fi Sco and 
fj, 2 Sco. 

Transformation s for Vt and Hp magnitudes to the 
V band are given bv lBesselll (|2000t) . |Platais et all (|2003l 1 also 
give a transformation of Vt to V as a function of B — V . Us- 
ing these two sources it was calculated that any correction 
for the colour difference between fj, 1 Sco and (i 2 Sco would 
be < 0.001 magnitude. Zero-point corrections of 0.041 and 
0.082 were thus applied to the Vt and Hp magnitudes re- 
spectively. 

The various systems used will have different internal 
errors, however a typical value for the errors of transformed 
photoelectric measurements is ±0.01 magnitude while errors 
for photographic photometry a re typically ±0.05 magnitude 
(|Budding and Demircanlll973l 1. 



4.2 Determining Effective Temperature from 
photometric indices or spectral types 

While separate spectra for the components of fi Sco can be 
discern e d, the ir spectral types remain somewhat uncertain. 
iMaurvl <| 19201 ) first assigned a spectral type of B3 to both 
components but noted a diff erence in their lin e inten sities 
which is discussed in detail bv lStruve and Elvevl [| 19421 The 
GCVS cl assifies the component s of u 1 Sco as B1.5V ± B6.5V 
although lStickland et alj (1996) discuss the possibility of the 
secondary being of a type earlier than B6. 

As the measured colour indices for fj, 1 Sco are for the 
combined light of the two components, tables of effective 
temperature and absolute magnitude as a function of spec- 
tral type or colour index cannot be readily used. It is, how- 
ever, useful to examine if suitable values for the temper- 
atures of the two components can be estimated from the 



combined-light colour index so as to assist with the initial 
choice of these parameters fo r subsequent modelling. 

Table |3] lists the GCPD j|MermiHiod et al.lll997r i mean 
photometric values in the UBV and uvby systems for /x 1 Sco 
(accessed online, 18 December 2008). Being at a distance 
of about 150 pc, interstellar reddening would be expected 
to be small for fx Sco. Using the Q value calculated from 
the UBV photometry, assuming a standard reddening law 
can be applied to the combined li ght of the components, 
and u sing the relationships given bv lBudding and Demircar] 
(|l973h gave an E B -v = 0.02. Additionally, from the uvby 
data, an E b _ v — 0.014 was calculated which is consistent with 
the calculated E B -v value. A value of (B — V) = —0.232 was 
thus adopted for the combined B — V colour index of fi 1 Sco. 

With the data in Table U (assembled from various 
sources), the B — V colour arising from combination of var- 
ious spectral types was then explored using a spreadsheet. 
Within the stated error of the listed B — V, the modelling 
indicated that the primary could not be earlier than B1.5V 
or later than B2V. For a primary of B1.5V the colour index 
varied only slightly as the secondary was varied in type from 
B3V to B8V. Within the listed error for the colour index, the 
spectral type of the secondary thus remains indeterminate. 
Best estimates for the effective temperatures of the compo- 
nents of u Sco using the combined B — V colour index, the 
spreadsheet modelling and the values given in Table [4] are 
thus 22,800 K for the primary and between about 12,000 K 
and 17,000 K for the secondary. 

The error in T etf may be estimated from the errors in 
log(T of f) for the stars used by iFlowerl (1 19961 ) to determine 
his B — V, T e tt relationship. For the 26 stars listed with 
-0.30 < B — V < -0.15 a mean error in log(T cff ) of 0.02 
was calculated. This translates to an error of approximately 
±1000 K for an early B star. 

Using the p a rallax of 6.51 ± 0.91 mas for /i 1 Sco given 
by Ivan Leeuwenl l|2007l ). along with a V ~ 3.0 for the com- 
bined light of the components, an M v ~ —2.9 ± 0.3 was 
calculated. This is consistent with a spectral type for the 
primary of B1.5V to B2V and in the range B3V to B8V for 
the secondary. 



4.3 Choice of Epoch 

The earliest published photoelectric measure- 
ments with an a ccomp anying epoch are those by 
iRudnick and Elvevl (| 19381 ) where they applied a cor- 
rection o f -0.0 06 day to Ma ury's published epoch . 

Ivan GentJ l| 19391 ). IStibbj l|l94ctl . IStickland et~ai1 j|l996ri . 
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Table 3. GCPD photometric values (mean) for fj, 1 
UBV and uvby systems. 



Sco in the 



V 

B-V 
U - B 
b~y 
m\ 
ci 



3.044 ± 0.023 

-0.212 ± 0.007 

-0.850 ± 0.005 

-0.089 

0.078 

0.103 



Table 4. Colours, absolute magnitudes and effective tempera- 
tures for main se quence B st ars. B — V indices a re fro m lZombecla 
lll99rj ), M v from lCoxl fcOOOll, T^t fro m lFloweil lll996ft , and b - y 
and M v (b - y) are from lMoor] l|l985h . 



Sp. Type 


B-V 


M v 


T eff 


b-y 


M v (b-y) 


B0 


-0.30 


-4.00 


33620 


-0.120 


-3.75 


Bl 


-0.26 




24338 


-0.115 


-3.31 


B2 


-0.24 


-2.45 


21261 


-0.092 


-2.08 


B3 


-0.20 




16958 


-0.080 


-1.33 


B5 


-0.16 


-1.20 


14203 


-0.068 


-0.97 


B6 


-0.14 




13197 


-0.059 


-0.74 


B7 


-0.12 




12368 






B8 


-0.09 


-0.25 


11376 


-0.045 


-0.32 


B9 


-0.06 




10612 


-0.037 


0.22 



iDanilekiewicz-Krosniak and Kurpinska-Winiarskal (|2002l ) 
and IShobbrookl (|2004 ~) also supply epochs based on the 
primary minimum. Table [5] summarizes the various epochs 
listed. Shobbrook's epoch was selected as a starting point 
as not only was it the most recent, but a small adjustment 
in the order of -0.001 day (within Shobbrook's stated 
error of ±0.002 day) resulted in a good alignment of the 
available datasets. By trying different values around that 
given by Shobbrook, an epoch for the primary minimum 
of HJD 2449534.17700(9) was selected as the value where 
there was the best alignment of the available photometric 
datasets. This was confirmed using the fitting routine 
within the PHOEBE software. 



4.5 The System's Orbital Characteristics 

PHOEBE was developed b y Prsa and Zwitter 



120031 ; IPrsa and Zwitterl 120051 ; IPrsa et al l |200S| ; IPrsal I200S 
and, while it is based on the Wilson-Devinney model 
(Wils on" and Devinnevl 197ll). it includes many re cent the- 



oretical developments ( Kallrath and Milond fl999l ) . Version 
0.31a of the model, (c. 2008), was used in this analysis. This 
software package requires separate radial ve locity data fo r 
each component; hence the extensive data of IMaurvl l|l920l ). 
which gives a radial velocity for the components combined, 
was not used. 

Using estimated spectral types for the components of 
H Sco, some parameters can be fixed in advance of mod- 
elling of the system's orbital characteristics. Owing to the 
radiative nature of B-type stars, the albedos of both the 
primary and secondary components were set to one as were 
the values for gravity darkening. The synchronicity param- 
eter for each star was also set to one. The option of an 
atmospheric model for the stars was enabled as was the op- 
tion of a reflection effect (set to four reflections). Factors 
such as third light, opacity, extinction and starspots were 
not considered. The limb darkening model ado pted is based 
on a logarithmic law and van Hamme's tables (|van Hammei 
I1993D . 

/j, Sc o is a member of the Scorpius- Centaurus OB As- 
sociation (jPreibisch and Mamaiek 2008 | ). th e nearest such 
association to the Sun. IP'Orazi et al.l 1 20091 ) have analyzed 
the metallicity distribution of open clusters within 500 pc of 
the Sun and note that the majority of the clusters with close- 
to-solar metallicity are part of the Local Association which 
includes the Scorpius-Centaurus association. The metallicity 
for the /x 1 Sco components was thus set to the solar value. 

The type of eclipsing binary to be modelled by PHOEBE 
is selected from a range of options. Initial attempts to fit 
the photometric data using some of the various options for 
the system failed to adequately represent features present 
within the data. It was found that the semi-detached op- 
tion, with the secondary filling its Roche Lobe appeared to 
provide a solution that best represents all the photometric 
data. The analysis thus concentrated on this option although 
a detached option was also investigated as it was able to rep- 
resent many of the key features as the size of the secondary 
approached its Roche limit. 



4.4 Determination of the Period 

The period l isted in the GCVS is based on that given by 
IStibbs] (ll948T). More r ecen t values are , howe ver, available 
from Shobbrookl l|2004l ) , or lArias et al l (|2005l ) . Table \5\ de- 
tails the various periods given in the literature with quoted 
errors where available. 

Starting with Shobbrook's period estimate, small ad- 
justments were made until the datasets aligned. It was 
found that a period of 1.4462700(5) days best represents the 
512 measurements collected over the period of the available 
measurements which is about 70 years. Figure [T] shows that 
there is excellent alignment between the various datasets 
using this refined value of the period along with the chosen 
epoch. The fitting routine within PHOEBE confirmed this 
result. 



4-5.1 Mass Ratio, Masses, Semi-Major Axis and 
Inclination 



Using the radial v elocity data of IStickland et al.| fl 996). 
based on IUE data ^International Ultraviolet Explorer Data! 
2009), along with the PHOEBE software, radial velocity 
curves were obtained and the mass ratio, semi-major axis 
and inclination determ ined. These val u es (sl ightly differ- 
ent to those found by IStickland et al.l (|l996t )L are given 
in Table [6] Figure [3] displays the radial velocity data from 
IStickland et alj (|l996l ) and resulting theoretical fit. 

While the values de rived here are consistent with those 
of IStickland et all 1 19961 ) to within the measurement errors, 
they are significantly different from those determined in ear- 
lier studies by Maury, van Gent, Rudnick and Elvey and 
Stibbs. 
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Table 5. Published epochs of primary minimum and periods. lArias et al.l (|2005h values based on 
iDanilekiewicz-Krosniak and Kur pinska-Winiarska (2002). 



Reference 



Epoch of Minimum (HJD) Period (days) 



Maury ( 19201 
Rudnick and Elvc 
van Gent (19391 
Stibbs (1948) 
Stickland et al. (1996) 
Shobbrook (2004) 
Arias et al. (2005) 
Proposed in this paper 



2412374.434 

2428281.250 

2428414.2978 

2432001.0451 

2448102.521 

2449534.178 

2432001.0475 

2449534.17700 



1.44627 
1.44627 

1.4462683±0.0000004 

1.44627±0.000001 

1.446271±0.0000044 

1.446270±0.000001 

1.44626876 

1.4462700±0.0000005 



Table 6. fi Sco system properties determined from modeling and analysis. 



Property 




Semi- Detached Model 


Primary Star's 


Spectral Type 


B1.5V [22,800 ± 1000 K] 


Primary Star's 


Mass 


8.49 ± 0.05 Mq 


Primary Star's 


Radius 


4.07 ± 0.05 R Q 


Primary Star's 


T cff (PHOEBE) 


23,725 ± 500 K 


Primary Star's 


Potential 


3.85 ± 0.01 



Secondary Star's Spectral Type 


B8 - B3 [12,000-17,000 K] 


Secondary Star's Mass 


5.33 ± 0.05 M 


Secondary Star's Radius 


4.38 ± 0.05 R 


Secondary Star's T eff (PHOEBE) 


16,850 ± 500 K 


Secondary Star's Potential 


3.07 ± 0.01 


Orbital Inclination 


65.4 ±1° 


Semi-Major Axis 


12.90 ± 0.04 R 


Eccentricity 


0.0 


Orbital Period 


1.4462700(5) day 


Centre of Mass Velocity 


-6.26 ± 0.04 km/s 


Mass Ratio 


0.627 ± 0.004 



Phase 

-0 60 -0.40 -0 20 00 20 




All Data 
-ITieoretical Frt 



Figure 3. PHOEBE's theoretical fit to all of the measured 
magnitudes of /i 1 Sco as function of phase using an epoch of 
HJD 2449534.1770 and a period of 1.4462700 day. 



4-5.2 Eccentricity, Argument of Periastron and Surface 
Potentials 

Using the available data noting the phase difference between 
the primary and the secondary minimums and the durations 



of the primary and secondary eclipses it is possible to de- 
termine the eccentricity. After normalizing the photometric 
datasets, a best-fit analysis using PHOEBE derived an ec- 
centricity of zero thus indicating the orbit is circular. Con- 
sequently the argument of periastron was set to zero along 
with the first time derivative of periastron. The circular or- 
bit is consistent with previous analyses of fi 1 Sco and that 
tidal interactio ns in close binaries tend to circulari ze the 
orbit over time (|HutJll98ll ; iKallrath and Milondfl999l) . 

The surface potentials were estimated using the mass- 
ratio, radii and semi-major axis values established from fit- 
ting the radial velocity curves. 



4-5.3 System Properties 

The (i 1 Sco system was modelled iteratively starting with the 
recent data acquired by M oon (presented in this paper) and 
that of IShobbrookl |2004) as they have the smallest photo- 
metric errors and together comprise half the the photomet- 
ric measurements taken to date. Using initial estimates pro- 
vided by this first phase, all available data was then used in 
the modelling of the fj} Sco binary system; the system prop- 
erties thus determined are given in Table [6] These values 
differ markedly from those determined by Maury, van Gent, 
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Figure 4. A geometric configuration for fj, 1 Sco at phase 0.25 
based on the modelling and analysis undertaken, axes in units of 
semi-major axis. 



axis, masses and radii, further data could be useful in im- 
proving the modelling of this binary system. Additionally, 
photometric measurements near phases 0.1, 0.4, 0.6 and 
0.9 could aid refinement of the current theoretical model 
through better defining some apparent features. 

With the semi-major axis of the orbit subtending 
0.39 mas, this star may be a suitable target for the Syd- 
ney University Stellar Interferometer (SUSI). 
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Rudnick and Elvey and Stibbs. In particular the masses cal- 
culated are significantly less. 

Similarly, the effective temperatures determined by 
PHOEBE were somewhat different to those determined in 
earlier studies but were found to be consistent with those es- 
timated from the combined colour index. The theoretically- 
derived fit, using the parameters in Table [6] and all available 
photometric measurements, are plotted in Figure [3] A geo- 
metric configuration for the stars at Phase 0.25 is presented 
in Figure [4] 

From the revise d parallax for this star of6.51±0.91 mas 
(van Leeuwen|[2007l) , the semi-major axis of the orbit is cal- 
culated to be 0.39 mas. This star may thus be a suitable tar- 
get for the Sydney University Stellar Interferometer (SUSI). 



5 CONCLUSIONS 

The light variations of /i 1 Sco appear to be stable over 
the 70 years for which photometric data exists. Using the 
comprehensive dataset of more than 500 photoelectric mea- 
surements assembled here, the period was determined to be 
1.4462700(5) day. A suitable epoch based on the primary 
minimum is suggested to be HJD 2449534.17700(9). 

Using all readily available data, and the PHOEBE soft- 
ware package, the data for fi 1 Sco could be best represented 
by considering the system to be a semi-detached binary 
in which the size of the secondary is close to, or fills, its 
Roche Lobe. The system properties determined from com- 
prehensive, iterative modelling using PHOEBE indicate that 
the masses of the components may be somewhat less than 
previously thought. With the secondary appearing to fill 
its Roche lobe it maybe possible to obtain spectroscopic 
evidence of gaseous streams in the /i Sco system noting 
that Doppler tomography has pr oduced indirec t images of 
gas flows in interacting binaries (Richards 2006). However, 
for Algol-type binaries the gaseous streams are faint rela- 
tive to the main sequence primary star and their detection 
challenges current techniques. Improvements to techniques 
would likely be required to detect such gas flows in /i Sco. 

There are currently only limited radial velocity data 
available for the fi 1 Sco. Given that the radial velocity data 
provides initial values for parameters such as the semi-major 
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